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We calculate the critical exponent 7 of Dyson's hierarchical model by direct fits of the zero 
momentum two-point function calculated with an Ising and a Landau-Ginzburg measure, and by 
linearization about the Koch-Wittwer fixed point. We find 7 = 1.299140730159 ± 10~ 12 . We 
extract three types of subleading corrections (in other words, a parametrization of the way the 
two-point function depends on the cutoff) from the fits and check the value of the first subleading 
exponent from the linearized procedure. We suggest that all the non-universal quantities entering 
the subleading corrections can be calculated systematically from the non-linear contributions about 
the fixed point and that this procedure would provide an alternative way to introduce the bare 
parameters in a field theory model. 



I. INTRODUCTION 



Scalar field theory has many important applications in 
condensed matter and particle physics. Unfortunately, 
there exists no approximate treatment of this theory 
which could pretend to compete in accuracy with pertur- 
bative methods in quantum electrodynamics at low en- 
ergy. Accurate calculations of subtle effects at accessible 
energies provide a window on hypothetical high energy 
degrees of freedom which are not accessible by produc- 
tion experiments. If we imagine for a moment that the 
kaons were the heaviest particles that we could produce, 
a precise determination of their weak matrix elements 
would become a unique way to obtain a quantitative in- 
formation about the charmed quark. Rescaled versions 
of this imaginary situation may become relevant in the 
future. 

The main goal of this article is to demonstrate that 
the use of hierarchical approximations allows determina- 
tions of the renormalized quantities, with an accuracy 
which can compete with perturbative QED, and for a 
wide range of UV cutoff and bare parameters. The use 
of hierarchical approximations simplifies the renormal- 
ization group (RG) transformation while preserving the 
qualitative features of scalar field theory. Well-known 
examples are the approximate recursion formula derived 
by K. Wilson Q, or the related recursion formula which 
holds for Dyson's hierarchical model gj. If used for quan- 
titative purposes, hierarchical approximations need to be 
improved. This is a difficult task under investigation 
and which is not discussed here. If the hierarchical ap- 
proximation could be improved in a way which maintain 
the advantages of the approximation, one could obtain 
results with an accuracy which would outperform any 
Monte Carlo method and defy any experimental patience. 

In the following, we use the scaling laws [|| to express 



the renormalized quantities as function of the bare quan- 
tities and a UV cu t-off. This parametrization (see e.g., 
Eqs. (|l.4| ) and (1.6)) is motivated below. We show with 
an example that the unknown parameters entering in the 
scaling laws can be determined accurately. In particu- 
lar, one universal quantity entering in the scaling law 
associated with the two point function (the critical expo- 
nent 7) can be calculated with two independent methods 
with an agreement in the 12-th decimal point. These 
two methods were sketched in Ref. Q. In the meantime, 
these methods were improved in order to get a signifi- 
cantly better accuracy. A detailed description of these 
two methods is the main technical content of the present 
paper. 

All the calculations reported in this paper were made 
with a specific example selected for its simplicity. How- 
ever, there is nothing essential in the choice of this exam- 
ple and accurate calculations can be performed with the 
same tools in a much broader range of models and param- 
eters. For the sake of definiteness we now describe the 
example chosen hereafter. We have limited the discussion 
to a calculation of the zero-momentum two point func- 
tion of Dyson's hierarchical model with a one-component 
scalar field and with a choice of the free parameter (de- 
noted c) which approximates a D = 3 theory. We consid- 
ered a wide range of UV cutoff (14 orders of magnitude) 
and two different sets of bare parameters. The choice 
of the hierarchical model is not essential cither. Wilson's 
approximate recursion formula is closely related to the re- 
cursion formula appearing in Dyson's hierarchical model 
[pj. It is possible to continuously interpolate between 
the critical exponents of the two cases [§. The numer- 
ical treatment of the two cases is completely identical. 
We have specialized the discussion to the case of Dyson's 
model because this model has been studied §-[[(]] in great 
detail in the past and because a fixed point of this model 
at D = 3 is known with great precision M. For the sake 
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of completeness the main features of Dyson's model are 
reviewed in section II A . The choice of the two-point 
function is not essential. The methods can be extended 
to other renormalized quantities as explained in Ref. [fl0| . 

The goal of a typical field theory calculation, is to 
obtain the renormalized quantities corresponding to the 
bare parameters entering in an action and a given UV 
cut-off A. In the following calculations, the bare pa- 
rameters will appear in a local measure of the Landau- 
Ginzburg (LG) form: 



(1.1) 



The UV cut-off corresponds to the scale where the theory 
under consideration stops being an accurate description 
and a more complete or more fundamental theory is re- 
quired. Given a set of bare parameters and a UV cutoff 
A, one can try to integrate M the degrees of freedom 
between A and a lower energy scale of reference Ar in 
order to obtain an effective theory describing phenomena 
at scales below A^. As an example, if we are interested 
in low energy processes involving pions, A could be cho- 
sen around m p and A^ around m n . This gives a ratio of 
about 6 between the two scales. Similar ratios may be 
applicable for an effective description of the Higgs boson 
in the hypothetical case that it results from an under- 
lying strongly interacting theory. If we are interested in 
the effects of the charmed quarks in non-leptonic decays 
of kaons, the ratio of the two scales would approximately 
be 3. Larger ratios appear if, for instance, we are inter- 
ested in the effects of the top quark in systems involving 
only the five other quarks, or the effects of the W and Z 
gauge bosons on the propagation of an an electron in a 
constant magnetic field. From these examples it is clear 
that one would like to be able to cover a wide range of 
values of Aji/A. 

In the scalar theory under consideration here, the cut- 
off is lowered by discrete steps which reduce the initial 
cutoff by a factor 2~"d for a "D-dimensional" theory (the 
noti on of dimensionality is explained at length in section 
|ll C| ). The limit of a large UV cut-off A can be reached 
by fine-tuning /3, the inverse temperature in Dyson's for- 
mulation of the model [|| . We use here the statistical me- 
chanics language: the magnetic susceptibility x is stud- 
ied by varying 0, keeping the bare parameters fixed. We 
could have, in a completely equivalent way, called the sus- 
ceptibility the zero momentum two point function, set 
- the kinetic term coupling constant - equal to one and 
fine-tuned the bare mass. 

We now follow Ref. jt) and consider a sequence L = 
1,2... of models with /3 = (/3 C — X[ L n) where Ai is the 
largest eigenvalue of the linearized renormalization group 
transformation and /i an arbitrary positive parameter. If 
we consider a fixed value of L and if \i is of order one, 
it takes about L iterations of the renormalization group 
transformation to get an effective theory with a mass of 
order one in A/j units. This come s fro m the fact that in 
the linearized approximation (see IIP ) , (3 C — (3 measures 



how far we are away from the stable manifold p] and 
this quantity is multiplyed by Ai at each iteration until 
it reaches a value of order one and the linearization does 
not hold any more. This suggests the definition of the 
renormalized mass m%: 
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where A^ is the UV cut-off defined as 



Ar = 2oA 



R ■ 



(1.2) 



(1.3) 



For (3 close enough to (3 C (i.e , for A large enough), one can 
approximate the susceptibility with an expression which, 
when D < 4, takes the form H 



x*{pc- PrHA, + Alfa - (3) A + . . .) 



(1.4) 



where Aq, Ax,.. . are functions of the bare parameters 
only. From t he ab ove equations and the expression of 7 
given in Eq. (2.17) which implies 



one obtains 



A 7 = 2° . 
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(1.5) 



(1.6) 



The main technical endeavor pursued in this article is 
t o d etermine numerically the unknown quantities in Eq. 
(1.4) and to determine the nature of the next corrections. 
We have used two independent methods. The first one 
consists in fitting the susceptibility at various values of 
/3 using Eq. (|1.4|). The second method consists in cal- 
culating the eigenvalues of the linearized renormalization 
group transformation in order to determine the critical 
exponents. The present article contains the details of the 
results announced in Ref. In the meantime we have 
refined some of the procedures used previously and im- 
proved significantly the accuracy of our results (e.g., at 
least four more significant digits in 7). These refinements 
are reported in the present article. 

The first estimation of the unknown quantities is based 
on a method of calculation presented in Ref. |l(| where it 
is shown that the use of polynomial approximations in the 
Fourier transform of the recursion relation allows efficient 
and highly accurate calculations of the zero-momentum 
Green's functions in the symmetric phase. The method 
is reviewed in section || where we also justify the in- 
troduction of a dimensionality parameter and review the 
linearization procedure. In section III, we analyze the er- 
rors associated with the method. We first explain how to 
get rid of the volume effects for the range of temperature 
considered later. We then analyze the round-off errors 
in arithmetic operations and show how to reduce them 
to an acceptable level by the use higher precision arith- 
metic, when necessary. We then discuss the effects of the 
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polynomial truncations and of the numerical errors in the 
calculati on o f the Fourier transform of the measure given 
by Eq. (1.1) and show that they can be reduced to a 
level where they will play no role in the discussion which 
follows. We conclude the section with an explanation of 
why a great numerical accuracy can only be achieved in 
the symmetric phase. 



In section [V, we fit the su scep tibility with the four 
parameters appearing in Eq. (1.4), neglecting the next 
subleading corrections indicated by the + . . .. After four 
successive refinements, all based on reproducible linear 
fits, we obtain values of 7 with a numerical stability up 
to the 13-th decimal point. In section |v| we analyze 
the next subleading corrections and show that neglect- 
ing them affects slightly the 12-th decimal point in 7. 
All the calculations in these two sections have been done 
with two different measures and gave compatible results 
for the universal quantities (7 and A). 

The most efficient way to calculate the critical expo- 
nent 7 and A consists in using the linearized renormal- 
ization group transformatio n fl near a fixed point. This 
is done explicitly in section VI. In ref. pL we gave con- 
vincing arguments indicating the uniqueness of the non- 
trivial fixed point for a large class of theories. In this 
article, we take this uniqueness for granted and we use 
the very accurate expression of this fixed point obtained 
from the work of Koch and Wittwer JjJ rather than the 
less accurate fixed points used in Ref. 0. This allows 
us to obtain values of 7 with estimated errors of less 
than one in the 13-th significant digit, the actual value 
agreeing with the previous estimate within the expected 
uncertainties. 

If there is only one non-trivial fixed point (universality) 
and if we can calculate accurately the exponents, the task 
of calculating the renormalized quantities for a particu- 
lar set of bare parameters reduces to the determination of 

quantities like Aq, A\, This task can be achieved by 

repeated subtractions as shown in sections IV and [v|. We 
are convinced that such calculations could be performed 
more efficiently by using the fixed point and a calculation 
of the nonlinear effects. We have checked [jllj that such a 
calculations can be satisfactorily performed in simplified 
versions of the basic recursion relation. If this task can 
be successfully completed for the model considered here, 
this would mean that the precise knowledge of the fixed 
point provided by Ref. is equivalent to a solution of 
the model. 

An accurate determination of the universal exponents 
provides a new approach of the renormalization proce- 
dure: we could try to treat as much as possible of the Ao, 
A\ , . . . and the corresponding quantities for the higher 
point functions as input parameters. This of course sup- 
pose that we have a detailed knowledge of their relative 
dependences. This question is under investigation with 
various methods. 



II. THE MODEL AND THE CALCULATION OF 
THE SUSCEPTIBILITY 

In this section, we describe the method used to calcu- 
late the magnetic susceptibility and introduce some def- 
initions which will be used later. For the sake of being 
self-contained, we first recall basic facts about Dyson's 
Hierarchical Model [0. 



A. Dyson's Model 

The model has 2™ max sites. We label the sites with 
n-max indices in m „,-,ii, each index being or 1. In 
order to understand this notation, one can divide the 
2" m " sites into two blocks, each containing 2™ m< " -1 sites. 
If x nmaj . — 0, the site is in the first box, if a;„ max = 1, 
the site is in the second box. Repeating this procedure n 
times (for the two boxes, their respective two sub-boxes , 
etc.), we obtain an unambiguous labeling for each of the 
sites. 

The non local part of the action of Dyson's Hierarchical 
model reads 



*=-i£(z> n £ (£ 



„...xi)) 



(2.1) 



The index n, referred to as the 'level of interaction' here- 
after, corresponds to the interaction of the total field in 
blocks of size 2™. The constant c is a free parameter 
which control the decay of the iterations with the size 
of the boxes and can be adjusted in order to mimic a 
ZJ-dimensional model. This point is discussed in more 
detail below. 

The field <p( Xn ,...,xi) is integrated over a local mea- 
sure which needs to be specified. In the following, we 
will work with the Ising measure, Wo(<p) = 5(c/> 2 — 1) and 
the Landau-Ginsburg measure of the for m g iven in Eq. 
(1.1). The hierarchical structure of Eq. (2.1), allows us 
to integrate the fields while keeping their sums in boxes 
with 2 sites. This can be expressed through the recursion 
relation 
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(2.2) 



where C„+i is a normalization factor which can be fixed 
at our convenience. 



B. Polynomial Truncations 



Introducing the Fourier representation 
dk 



Jk<f> 



2tt 



W n (k) 



(2.3) 
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and a rescaling of the source by a factor 1/s at each it- 
eration, through the redefinition 



(2.6) 



Rn{k) = W n {\), (2.4) 
s 

the recursion relation becomes 

Rn+i(k) = C n+1 exp{~\p{-^T +1 ^){R n ^)f. (2.5) 

The rescaling operation commutes with iterative integra- 
tions and the rescaling factor s can be fixed at our con- 
venience. 

We will fix the normalization constant C n is such way 
that Rn(0) = 1. Then, R n (k) has a direct probabilis- 
tic interpretation. If we call M n the total field <fi x 
inside blocks of side 2 n and < ... >„ the average calcu- 
lated without taking into account the interactions of level 
strictly larger than n, we can write 



We see that the Fourier transform of the local measure 
after n iterations generates the zero-momentum Green's 
functions calculated with 2 n sites and can thus be used 
to calculate the renormalized mass and coupling constant 
at zero momentum. 

In the following, we use finite dimensional approxima- 
tions of degree l m ax of the form: 

Rn(k) = 1 + a rhl k 2 + a rh2 k 4 + ... + a n , lmax k 2lma * . (2.7) 

After each iteration, non-zero coefficients of higher order 
(a„-|_i.i maa . +1 etc. ) are obtained, but not taken into ac- 
count in the next iteration. More explicitly, the recursion 
formula for the a n ^ m reads : 



S;"m(S P+ ^°n, P an, g )[(20!/q-m)!(2m)!](c/4)'[-(l/2)^-"> 
Eb*(£ P+g=i a„, p a„, 9 )[(20!//!](c/4)'[-(l/2)/3]' 



(2.8) 



As one can see that once an initial Ro(k) is given, the 
procedure is purely algebraic. The initial conditions for 
the Ising measure is Ro{k) — cos(k). For the LG mea- 
sure, the coefficients in the fc-expansion need to be eval- 
uated numerically. This method has been discussed and 
tested at length in Ref. |10| . The dimension l ma x of the 
polynomial spaces required to make reasonably accurate 
calculation is remarkably small: less than 50 for a typical 
calculation (see Ref. [jl(| for details). 

As far as numerical calculations are concerned, the 
choice of s is a matter of convenience. For the calcu- 
lations in the high temperature phase (symmetric phase) 
not too close to the critical points, or for high tempera- 
ture expansions the choice s — \[2 works well H|[. On 
the other hand, the choice of rescaling factor s = 2c~ 1 / 2 
prevents the appearance of very large numbers when we 
are very close to the critical temperature. In the follow- 
ing, the finite volume magnetic susceptibility is defined 
as 



Xn(f3) 



< (M„) 2 > r 



(2.9) 



From Eq. (2.6), we obtain 



action given in Eq.(2.1) is invariant under the removal 
of the I = 1 terms (first level interactions) followed by 
the rescaling of the fields. In other words, the kinetic 
term is not renormalized and r/ =0. From this, we can 
derive the way c should be tuned in order to mimic a 
D-dimcnsional system. Given that the dimension of a 



scalar field in /^-dimension is [</>] = [L\ 
a length, we obtain in the continuum 



(£>-2) 



[( / d D xm) 2 \ = l 



D+2 



On the lattice this becomes 



[< (M n y >„] = l 



_ T D+2 



where L is 



(2.11) 



(2.12) 



If we use the rescaling factor s — 2c 1 / 2 , the non-local 



part of the action given in Eq. (2.1) is invariant under a 
renormalization group transformation. If in addition the 
local measure is also left invariant, the average values of 
the even powers of the rescaled field stays constant. Re- 
turning to the original field variables, we found that at 
(or sufficiently close to) a fixed point, 



-2an,i(y) r 



(2.10) 



< {M n y >„oc (-)" 

c 



(2.13) 



C. Introducing the Dimensionality 

From a conceptual point of view, the choice s = 2c -1 / 2 
is of particular significance because the infinite volume 



The only relevant scale is the size of the box over which 
we have integrated all the field variables except for their 
sum. The volume of the box is proportional to the num- 
ber of sites inside the box: 



L D oc 2 r 



(2.14) 



4 



Using this together with Eqs. (2.12) and (2.13) we obtain 



A. Volume Effects 



Z = 2 T>( D + 2 ) 
c 



(2.15) 



or in other words, c = 2 1 o . 

All the calculations done hereafter have been done for 
D = 3. 



D. Review of the Linearization Procedure 

We now briefly review the linearization procedure. We 
denote the eigenvalues of the linearized RG transforma- 
tion by A„ with the convention Ai > 1 > A2 > A3 . . .. The 
closeness to the fixed point is essentially monitored by the 
motion along the unstable direction. Until the number of 
iterations n reaches a value n* such that A™ (f3 c — [3) ~ 1, 
R n is "close" to the fixed point and a nj i stays close to its 
fixed point value (assuming that we use the scaling factor 
s = A=). When n gets larger than n*, x starts stabiliz- 



ing. Using the relation between a n ,i and Xn given by Eq. 
(|2.10p, we obtain the order of magnitude estimate: 



X -(-)"* =(/3c-/3) _ ™ 
c 



(2.16) 



Reexpressing in terms of {(3 C — (3) , we find that the expo- 
nent for the leading singularity is 



7 =ln(|)/ln(Ax) 



(2.17) 



According to the same linear argument, the order of mag- 
nitude of the components in the stable directions should 
be proportional to A™ with I > 2. Using the estimate for 
n* and reexpressing in terms of (0 C — (3), we obtain the 
subleading exponents A/ = — ln(A/)/ln(Ai) for I > 2. In 
the following we simply use the notation A for A2 and 
the higher exponents will play no significant roles. 



III. ERROR ANALYSIS 



There are three important sources of errors which need 
to be considered when we calculate the magnetic suscep- 
tibility: the finite volume effects, the round off errors and 
the effects of the finite dimensional truncation. A gen- 
eral discussion of these questions is given in Ref. Jh}. In 
the following, we discuss them in the particular cases re- 
quired for the calculations of section IV. In addition, we 



discuss the effects of the errors on the initial coefficients. 
All the calculations done hereafter have been made in the 
symmetric phase. In the last subsection, we explain why 
the present methods do not yield accurate results in the 
broken symmetry phase. 



As explained in II D| , when calculating the suscepti- 
bility at values of /3 close to and below (3 C , we spend 
about —ln(p c — P)/ln(Xi) iterations near the fixed point. 
During the se it erations, we have the "conformal" scal- 
ing of Eq.( 2.13 ) and the round-off errors are amplified 
along the unstable direction (see next subsection) . After 
that, assuming we are in the symmetric phase, the or- 
der of magnitude of the susceptibility stabilizes and the 
corrections get smaller by a factor | at each iterations. 
At some point, all the recorded digits stabilize (irrespec- 
tively of the numerical errors which occurred in the first 
stage described above). This give the estimate |l(J for 
the number of iterations n((3, P) to stablize P digits (in 
decimal notations) 



n(f3,P) = 



Dln(10) 
2ln(2) 



[P - jlog 10 (l3 c - 0)} . (3.1) 



For P = 16, 7 ~ 1.3 and (3 C — {3 = 1CP 9 , we obtain 
n ~ 140 and we need to add about 7 iterations each time 
we get closer to (3 C by a factor 1CU 1 . It is thus quite easy 
to get rid of the volume effects. In the following, we will 
perform calculation (3 < [3 C — 1CU 14 and n max = 180 will 
be enough to avoid finite volume effects. 



B. Numerical Errors 



From Eq.(2.S), we see that the calculation of each of 
the a n+ u involves a number of arithmetical operations 
proportional to l m ax- When we are close to the fixed 
point, these errors generate small contributions in the 
unstable direction. These errors are then amplified by a 
factor Ai at each iteration until we move sufficiently far 
away from the fixed point. Consequently, the closer (3 is 
to (3 C , the more time is spent near the fixed point and 
the larger the numerical error become. A simple calcula- 
tion |10) corresponding to this reasoning shows that the 
relative errors obey the approximate law 



(3.2) 



where S is a typical round-off error. 

A simple way to probe the numerical errors is to make 
a small change in the rescaling factor s. As explained in 
the previous section, we can in principle use any value 
of s to calculate the susceptibility. This arbitrariness is 
compen sated at the end by an appropriate rescaling given 
in Eq. ( 2.10 ). If we could perform the arithmetic opera- 
tions exactly, the susceptibility would be completely in- 
dependent of s. However, due to the round-off errors, the 
susceptibility actually depends on s. This is illustrated 
in Fig. [l] where we calculated the distribution of x f° r 
values of s varying between 

by step of 10~' 
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FIG. 1. Distribution of the magnetic susceptibility x with 
respect to the scaling factor s. 

This calculation has been performed in Fortran with 

double precision variables. We have used an initial Ising 

measure with (3 = (3 C — 10~ 9 . This distribution has 

a mean value /i = 1.041926904 x 10 12 and a variance 

a = 4.9 x 10 5 . From these quantities, we estimate that 

the relative errors I — I due to numerical errors should be 
1 x 1 

of order <jj \l = 4.7 x 10~ 7 . This is in agreement with the 



order of magnitude estimate of Eq. (3.2): using 10 as 



a typical round-off error in a double-precision calculation, 
we obtain |-^| ~ 10~ 7 for (3 C — (3 = 10~ 9 . A more accu- 
rate calculation performed with methods described below 
gives the result x = 1-041926626 x 10 12 . We have checked 
that this result was invariant under slight changes in s. 
From this, we see that /i — \ — 2-8 x 10 5 which is ap- 
proximately 0.57a. Another information concerning the 
spread of the values is the difference between the largest 
and the smallest values of the distribution which is 6.4a 
in the present case. 

There a several detailed features of this distribution 
which are not well understood. The first one is that the 
distribution is not symmetric about s — The values 



of x spread more above s = -4^. In addition, a more 
detailed study shows that the distribution is not well 
centered and that of values about the mean value de- 
parts more from a gaussian distribution than expected, 
given the number of "independent trials" (2000 in Fig. 
[TJ) made. In addition, increasing the statistics does not 
decrease n — x or increase significantly the difference be- 
tween the largest and the smallest value. These questions 
are now being investigated with low l ma x examples. 

In conclusion, we have a good control on the maxi- 
mal errors made as a consequence of the round-off er- 
rors. These seem not to exce ed 10 times the order of 
magnitude given by Eq. ( |3.2[ ). On the other hand, we 
have an incomplete understanding of their distribution 
within these bounds. These precludes the use of statis- 
tical methods to obtain more accurate results and other 
methods need to be used. 

The most efficient way to improve the accuracy of x 



consists in using higher precision arithmetic. This can 
be done easily, for instance, using the Mathematica en- 
vironment where one use the instruction SetPrecision[ ] 
to introduce numbers with a desired number of signifi- 
cant digits and the instruction Precision[ ] to monitor 
the numerical errors. The initial precision can then be 
adjusted empirically in order to obtain a desired accu- 
racy for x- This accuracy is then checked by making 
changes in s as explained above. A typical calculation 
with l max = 50, n m ax = 200 and a required accuracy of 
16 digits in the final result takes of the order of 10 3 sec 
on a common workstation. The same type of calculation 
in double-precision Fortran takes about 0.1 sec. While 
the high-precision program runs, we could thus run the 
double-precision 10 4 times. If a proper understanding of 
the statistical distribution of the errors was at hand (as 
explained above, this is not the case), we could hope to 
use the 10 4 values to reduce \5x\ by a factor 10~ 2 . In 
the example discussed above, we would get hope to get 
errors in the 9-th significant digit instead of the 7-th. 
However, with the high-precision method we obtained 16 
correct significant digits. In the example discussed above 
we obtain the accurate value x — 1-0419266255... x 10 12 . 
The difference between this more accurate value and the 
mean calculated above with 2000 data points is 0.57 x a 
and stays at the same large value when the statistics is 
increased. In the following, the high-precision method 
will exclusively be used. 



C. Determination of l„ 



As shown in Ref. IjlfJ, the effect of the finite trunca- 
tion decays faster than exponentially with l ma x in the 
symmetric phase. In general, the determination of l m ax 
depends on how far we are from criticality and the re- 
quired accuracy on the value calculated. In the following, 
we will require a 13 significant digits on x an d (Pc ~ P) 
with (3 < I3 C - 10~ 14 . When (3 = {3 C - 10~ 14 , all the 
significant digits up to the 13-th decimal point of the 
quantity (f3 c — (3) are lost since they cancel. Conse- 
quently, in order to get 13 significant digit in {(3 C — (3) 
in the range considered, we will determine (3 C with an 
accuracy of 10 -27 . This will be the most stringent re- 
quirement to determine l ma x- As explained in the previ- 
ous subsection, we can easily perform calculations with 
high-precision arithmetic and follow the bifurcations |Tc| ] 
in the ratios of successive a Ut i in order to determine (3, 
Fig. U shows the effect of adding or subtracting 10 
to (3 C = 1.179030170446269732511874097 in the case of 
an initial Ising measure. This calculation has been per- 
formed with l max = 50. If we use larger values of l m ax, 
(3 C remains at the quoted value. This is just a particular 
example. 

In general, the minimal value of l max for which (3 C 
stabilizes can be obtained from extrapolation from the 
changes at low l m ax where calculations take little time. 



27 



6 



We use the notation 5f3 c for (3 c {l m ax) — flc- The quan- 
tity logio|«5/3 c / /3 C | versus l ma x is shown in Fig. p| for the 
Ising model and the LG measure of Eq. (1.1) with 



= 1, p = 2, and g = 0.1. 
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A/3 C as a function of the relative 
-^f- = 10~ 2 (empty circles), 



FIG. 3. The shift in f3 c 
errors in the Z-th coefficient. = 10" 

10 -4 (filled boxes), 10 -6 (empty triangles) and so on until 

Sa O,l _ 1Q-14 
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FIG. 2. logio\5/3 c / (3c\ versus l ma x for the Ising case (cir- 
cles) and the LG case (stars). The solid line is a fit with 

a + b(lmaxln(lmax)) ■ 



The logarithm of the relative errors falls faster than 
linearly. In good approximation E3|, logio(A/3 c ) ~ 
a + b(l max ln(l max )). So if we want, say 5/3/ P ~ 10~ 27 , 
the choices of l max = 48 for the Ising case and l ma x = 60 
the LG cases appear to be a safe. In order to check the 
stability of these values, we increased the value of l ma x 
to 55 in the Ising case and to 64 in the (LG) case, and 
we obtained the same /3 C value in both cases. We have 



The results can be read as follows. If we are interested 
in determining f3 c with, say, 10 significant digits, 00,10 
has to be determined with 2 significant digits, ao,g with 
3 significants digits, ao.s with 4 significant digits etc... . 
In the following, we are interested in universal properties 
(features which are independent of the measure) rather 
than in properties of particular measures. Consequently 
we have only used a double-precision calculation of the 
Fourier transform for the LG model. The reproducibility 
of the details of the calculations then require having the 
same olqj. On the other hand, in the Ising case, the ana- 
lytical form of the initial coefficients allows a completely 
reproducible procedure. 



obtain 13 significant digit for \ in the range of (3 specified 
above. 



E. What Happens in the Broken Symmetry Phase 

Fig. U shows the existence of two phases. There arc 
five parts of the graph we would like to discuss here. 



D. Effects of the Errors on the Initial Coefficients 
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In the Ising case, Ro(k) = Cos(k), the initial coeffi- 
cients are known analytically: a$ : i — (— l) l /(2l)\. How- 
ever, this is not the case in general. We want to study 
the effect of a change in <5ao,; in the initial coefficients on 
f3 c . The results are shown in Fig. E| 
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FIG. 4. a n +i,i/ a ™,i versus n for f3 — /3 C =F 10 27 . 

For the rest of the discussion, it is important to specify 
that a n ^i has been calculated with the canonical value of 
the rescaling parameter s = The first part is before 
the bifurcation. This is shown as the region "1" in Fig. |] 
where the ratios of CL n +i,i/ a n,i are close to 1. The second 
part ( "2" ) shows the bifurcation in the high-temperature 
phase. If we are below f3 c , the ratio a n +i,i/^n,i w iU g° 
down to the value c/2, which guarantees the existence 
of a thermodynamical limit for \ (since we n eed t o mul- 
tiply a n ,i by (~)™ in order to get \ see Eq.( 2.1C )). On 
the other hand, the bifurcation toward the low temper- 
ature phase is characterized by a peak shown by "3" in 
Fig. [|. Part 4 of the graph is a narrow "shoulder". In 
the low-temperature phase, we expect < M 2 > rl cx 2 2 ", 
which means an+i^/a-n,! — c — 1-26. We studied the 
Imax dependence of this shoulder and observed that the 
number of points on the shoulder increases by approxi- 
mately one when we increase l max by 10. Unfortunately, 
the shoulder is not infinite and after a few iterations, the 
ratios will reach to 1 again (part "5"). This signals an 
attractive fixed point. However, this is not a fixed point 
of the exact (not truncated) recursion relation. This can 
be seen by looking at the coefficients a* of these attrac- 
tive fixed points for different values of l max - When l max 
increase, the values of a* increase like (l ma x) ■> showing 
that their existence is due to the truncation process. 

The fact that the truncation procedure generates nu- 
merical instabilities in the low-temperature pha se can be 
understood from the basic formula Eq. ( |2.5| ). In the 
low temperature phase, the measure W(cj>) has two peaks 
symmetrically located with respect to the origin. At each 
iteration, the separation between the peaks increases by 
a factor 2 (in unrescaled units). By taking simple ex- 
amples and going to Fourier transform, one sees that at 
some point the partial sums (t runc ated at l ma x) repre- 
senting the exponential in Eq. ( |2.5| ) becomes inaccurate 
because the argument of the exponential is too large. 



IV. CRITICAL EXPONENTS FROM FITS 



In this section, we explain how to calculate the expo- 
nents 7 and A using a sequence of increasingly accurate 
fits of the susceptibility. The general method has been 
briefly outlined in Ref. Here, we give all the de- 

tails of a significantly more accurate calculation which 
leads to a determination of 7 with 12 decimal points. 
The main ingredient of the procedure is that for (3 close 
enough to (3 C (i.e, for a cut-off A large enough), one can 
approximate very well the magnetic susceptibility (zero- 
momentum two point function) with an expression taking 
into account only the first irrelevant direction, namely : 



The estimation of the unknown quantities in this equa- 
tion proceeds in four steps. In the first step, we get a 
rough estimate for 7 by using a linear fit in a range of /3 
where we minimize the combined effects of the numerical 
errors and of the subleading corrections. In the second 
step, we discuss how to improve this result by estimating 
the sub-leading exponent A and the coefficient Ai/Aq. 
Using these preliminary estimates, we will as the third 
step of the procedure, use a "bootstrap" technique be- 
tween a set of high-precision data close to criticality and 
another set of data where the subleading corrections are 
important. Finally, we do a linear analysis of the differ- 
ence between the fit and the high-precision data in order 
to get results which are as independent as possible of the 
slightly arbitrary choices (how to divide the data into 
"bins" etc... ) made during the first three steps. After 
this fourth step, we analyze the difference between the fit 
and the data away from criticality and discuss the next 
subleading corrections. 

All the calculation of these sections have been done 
with cither an Ising init ial m easure or a LG measure of 
the form given in Eq. (1.1) with m 2 = 1, p = 2, and 
g = 0.1. We refer to these choices as the "Ising case" or 
the "LG case" hereafter. 



A. Localized Linear Fits 



In the first step, we calculate x at various tempera- 
tures and display logio(x) versus -logio(/3 c — /?)■ We will 
use the notation 



x = -log w ((3 c ~ P) 



(4.2) 



If we display logio(x) versus x , we see a linear behavior 
with a slope 7 ~ 1.30. The deviations from the lin- 
ear behavior are not visible to the naked eye. We need 
to study these deviations locally in /3. In order to un- 
derstand the corrections, we have divided the data into 
14 bins of 100 points. The first bin contains the data 
x = 1.00, 1.01, 1.99 and so on. In each bin, indexed i , 
we make a linear fit of logio(x) versus x. In the i-th bin 
we will call the slope and (er^) 2 denoted the sum of 
the squares of difference between the data and the lin- 
ear fit divided by the number of points in a bin minus 2 
which is ( for the i-th bin ), 



(a 



(*h2 _ 



data\ 



io 9w {xi:-)f 



()<s 



(4.3) 



X~ (/3 C -/3)^(A + A 1 (/3 C -/5) A ) . 



(4.1) 



where j indexes the data points in the i-th bin. The 
values of logio^W) are plotted in Fig. |[ 
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FIG. 5. The deviations from the linear fits logioo"' 1 ' defined 
in the text as function of the bin index i, for the Ising model 
(circles) and the Landau- Ginsburg model (stars). 

It is easy to interpret this graph. There are two ma- 
jor sources of deviations from the linear behavior. The 
first one is the existence of subleading corrections to the 
scaling laws which decrease when j3 gets closer to f3 c . As 
a first guess, we use cr") cx 10 _Aj so that logioc^' — 
— Aj + constant. By calculating the slopes between bin 
1 and bin 9, we obtain A ~ 0.42 and A ~ 0.45 for the 
Ising and the LG cases, respectively. Thus, we already 
obtained a numerical value for the subleading exponent 
which is roughly the same for the two models considered 
here. The other source of deviation from the linear be- 
havior comes from the numerical errors discussed in the 
previous section and which increase when [3 gets closer 
to (3 C according to Eq. (3.2). a ~ 6/((3 c — 0). The slopes 
between bin 9 and bin 14 are —0.95 and —0.96 for the 
Ising and the LG cases, respectively, in good agreement 
with the (p c — P)^ 1 dependence of the numerical errors 
predicted by Eq. (fTs]). 

In bin 9, these two deviations from linear behavior are 
minimized and we can consider 7^) as a first estimate 
of 7. Its numerical values is 1.29917 for the Ising case 
and 1.29914 for the LG case. By using this simple pro- 
cedure we already gained almost two significant digits 
compared to the existing estimated |^,^| where the an- 
swer 7 = 1.300 was obtained with errors of order 1 in the 
last digit. 



B. Subleading Corrections 

The second step consist in correcting the previous esti- 
mate by taking into account the subleading corrections. 
We will use the bins 6 and 7 where the next subleading 
corrections are reasonably small and the numerical errors 
are not too large. We have divided these two bins into 
10 sub-bins of 100 points. We will use two digit indices 
for these sub-bins. For instance, sub-bin 6.5 is the fifth 
sub-bin of bin 6 and contains the values of x: 6.5, 6.501, 



... ,6.599. Using the notation x for the middle of the 
sub-bin and Eq. (LI) we obtain, 



log w {x{x)) =ix + log w (A + AilO Aa; ) 



(4.4) 



For a small change Sx with respect to x, we obtain that 
at first order in this change 



logio{x{x + Sx)) - log 10 (x(x)) 

(7- ^AW- As )5x + ${6x 2 ) . 
Aq 



(4.5) 



The coefficient of 5x can be interpreted as the local slope 
near x. Indexing each sub-bin by j (e.g j : 6,6.1,6.2...) 
and its middle by j + 0.0495 (e.g. 6.1495 is the middle of 
the sub-bin 6.1), the the slope yV' in the sub-bin j reads 



7 (j) ~ 7 _ A(— )1 - A «+0.0495) 
Ao 



(4.6) 



The unknown quantities ^ and A can be obtained 

from linear fits of logio(|7^ +0 ' 1 ^ — 7 I) versus j. From 
Fig. [| one can see that in good approximation, there ex- 
ists an approximate relationship between these two quan- 
tities in the two cases considered. In addition, the slope 
appears to be identical for the two cases. 
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FIG. 6. The linear fits of logi (|7 0+0 - 7 (j) |), for the 
Ising case (circles) and for the LG case (stars). 



Th e slope and the intercept can be calculated from Eq. 
(4.6) which implies that 



log w {\l U+0 - 1} ~ 1 U) I) - -A(j + 0.0495) + 

z 03lo (i-io-°- 1A ) 



^0 



(4.7) 



The subleading exponent A is the absolute value of the 
slope. Having determined A and knowing the intercept 
we can then determine Zo<7io(|j|j|). Using this procedure, 
we obtained Ai/A = -0.57 and A = 0.428 for the Ising 
case and Ai/A Q = 0.14 and A = 0.427 for the LG case. 



9 



If we now repeat the first step - a linear fit in bin 9 
- but with x divided by [1 + (Ai/A Q )(l3 c - /3) A ], we ob- 
tain 7 = 1.299141 =F 10~ 6 for the two models considered 
above. Given that for a calculation using double preci- 
sion in bin 9 we have values of x with between 6 and 7 
significant digits, this estimate seems to be the best re- 
sult we can obtain with this procedure (see the discussion 
of the numerical errors in section III.B). 



C. A Bootstrap Procedure Involving Higher 
Precision Data 



Up to now, our desire to minimize the subleading 
corrections which decay like 10~ Ax has been contra- 
dicted by the appearance of numerical errors growing like 
10 (-i6+x)_ However, we have explained in section III.B, 
that it is possible to circumvent this difficulty by using an 
arithmetic having a better precision than the usual dou- 
ble precision. In this subsection, we will use data having 
at least 13 correct significant digits in bin 11, 12 and 13. 
We call this data "high precision data" . As we explained 
before, we choose l m ax = 50 and l ma x = 58 for the Ising 
and the LG case, respectively. 

Since the calculations are more lengthy, we used only 
10 points in each bin. We also determined (3 C with 27 sig- 
nificant digit so that in bin 13, the subtracted quantity 
(f3 c — (3) is also know at least with 13 significant digit. 
We found 0c = 1.179030170446269732511874097 for the 
Ising case and j3 c = 1.14352915687979895500964720 for 
the LG case. We then use bin 13 (where the subleading 
corrections are small and the errors are not very impor- 
tant) to calculate x divided by the subleading correc- 
tion as explained in the previous subsection to estimate 
7. Then with the new value of 7 obtained we go back 
to bin 7 to calculate the subleading corrections. This 
procedure can be iterated and this 'bootstrap' of lin- 
ear fits converges rapidly. We obtain 7 = 1.299140732, 
A = 0.4262 and A x /A = -0.564 for the Ising model and 
7 = 1.299140730,A = 0.4258 and A 1 /A = 0.135 for the 
LG case. These numbers change typically by one in the 
last digit quoted above if one replaces bin 7 by bin 6 to 
evaluate A. In order to remove this arbitrariness, we will 
now use these numbers as the initial values for a more 
accurate procedure. 



D. Linear Analysis of the Discrepancies 

We have now reduced the errors made in the estimate 
of the unknown quantities appearing in Eq. ( ]4. L| ) to a 
sufficiently low level to allow us to treat these errors in 
a linear approximation. We start with an initial fit of 
the data, for instance as obtained in step 3, with errors 
in the unknown quantities parametrized in the following 
way: 



[logiox]fit = logioA + S(log 1Q A ) + (7 + Sj)x + 

Z 0ffl o(l + ( C i+5 Cl )10-( A + 5A ^) • (4.8) 

where c\ = Ai/Aq and 7, A stand for the exact values. 
On the other hand, we assume that the data can be fitted 
according to Eq. (|4.1|) 



[logwxjdata = {log w A Q + jx + log w {\ + cilO x )) . 

(4.9) 

Combining the two above equations we obtain at first 
order 

[logwxlfit - [logwxUata ^ S(log 10 A ) + 5~fX + 



Zrc(10) 



- Cl 10- A;r x<SA . (4.10) 



Interestingly, the x-dependences of the four terms are all 
distinct and we can fit t>(logio^4o)> & c i an d SA us- 
ing a standard least square procedure where the function 
to be fitted depends linearly on the fitted parameters. 
This procedure can be repeated until some numerical 
stability is achieved. The final results are insensitive to 
small changes in the initial values coming from the un- 
certainties associated with the previous step. Using bin 
13, we obtain 7 = 1.2991407301599 and A = 0.4259492 
for the Ising model, and, 7 = 1.2991407301582 and 
A = 0.4259478 for the LG case. The small numerical 
fluctuations which persist after many iterations produce 
changes of less than 2 in the last quoted digit. The ori- 
gin of these small fluctuations can be inferred by plot- 
ting Uogxox]fit - [logioxldata for the final fit (see Figs. @ 
and H in the next section). The non-smoothness of these 
differences in bin 13 indicates that they are due to the 
numerical errors on x- The amplitude of these differences 
is smaller than 10~ 13 consistently with the fact that we 
performed the calculation of % in a way that guaranteed 
at least 13 accurate significant digits. These fluctuations 
are indicative of the limitation in the numerical precision 
of our procedure. The accuracy of the value of the ex- 
ponents, i.e. how close they are to the "true" values, is 
further limited by the fact that there exist corrections to 
our main assumption Eq. (fl.l| ). If we assume universal- 
ity, the discrepancy between the values of the exponents 
for the two cases considered should give us an indication 
concerning the accuracy of the results. For instance, the 
discrepancy between the two estimates of 7 is of order 
10~ 12 which is about ten times larger than the fluctu- 
ations of numerical origin. The estimation of the next 
sub-leading corrections is the main topic of the next sec- 
tion. 



V. THE NEXT SUB-LEADING CORRECTIONS 

In the previ ous section, we have used the parametriza- 
tion of Eq. (4.1) for the susceptibility near criticality. 
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This parameterization is by no means exact and correc- 
tions become more sizable as we move away from criti- 
cality. The corrections come from effects which can be 
calculated by a linearization p roced ure (the next irrele- 
vant directions, see subsection II D) or effects which are 
intrinsically non- linear. Anticipating the results which 
will be presented in the next section, we obtain the ap- 
proximate values A2 ~ 0.43 (we recall that since we only 
took into account one irrelevant direction, we used the 
notation A for A2 before) and A3 ~ 2.1. In other words 
if we consider the first irrelevant direction as "first or- 
der" , the next irrelevant directions produce effects which 
are smaller than the fourth order. In bin 13, these effects 
are completely unnoticeable in our analysis. The non- 
linear effects are discussed in the subsection VA. The 



main result obtained there is that all the corrections can 
be parametrized in the following way: 

X ~ (& - /3r 7 (A + Ai((3 c - (3) A + A 2 (f3 c ~ /3) 2A 
+ A*{p c -0) + ...), (5.1) 

In the next subsection, we analyze the data in terms of 
the new parametrization and extrapolate our results in 
order to estimate the errors made in the calculations of 
7 and A in the previous section. 



we obtain the claimed 14 orders of magnitude. Such a 
small effect is smaller than our numerical resolution. 

Second, the singularity [fi c — /3)~ 7 should be replaced 
by {{fie ~ 0) + d 2 (P c ~ P) 2 + ■ ■ -)~ 7 with coefficients di 
calculable in low dimensional maps. These corrections 
generate analytical corrections to the scaling law in con- 
trast to the subleading corrections which are in general 
not integer powers 

Third the nonlinear corrections associated with the ir- 
relevant directions generate corrections which are pre- 
sumably of the form (/3 C — (3) lA with I = 2,3,.... Later 
we call these corrections the quadratic corrections or the 
second order effects. 

In summary, the corrections associated with n onli near 
contributions obey the parametrization of Eq. ( |5.l| ) for 
a sequence of exponents 0.43, 0.86, 1, 1.29, 1.72, 2,.... 
Note that these exponents are very close to each other 
and it may be difficult to disentangle their effects. 



B. Empirical Determination of the Corrections 



A. Nonlinear Corrections 

The previous analysis, describes the linearized flows 
near the fixed point. The closer to criticality we are, 
the more iterations are spent close to the fixed point and 
the more accurate the linear description is. Neverthe- 
less, when we approach or leave the fixed point, nonlin- 
ear effects are unavoidable. These nonlinear effects can 
be studied more easily in low-dimensional maps. With- 
out entering into the detail of this analysis [ pT| , we can 
envision three types of corrections which we now proceed 
to discuss. 

As already noticed in Q, the "constants" Aq and A\, 
should be replaced by functions Ai((j3 c — (3) such that 
Ai(Xi(f3 c — 0)) = Ai((3 c — (3). This invariance implies 
an expansion in integral powers of (/3 C — /3) lw (Fourier 
modes) with to = jjf-y ~ 17.8. We argue here that the 
coefficients of the non-zero powers are suppressed by 14 
orders of magnitude. The rationale for this suppression 
is that the non-zero modes contribute to the extrapo- 
lated slope (an asymptotic estimator for 7 — 1 used with 
the high-temperature expansion) with about the same 
strength as the zero mode ||. However, this is the re- 
sult of a double amplification for the non-zero modes. 
This results from the equations (3.7) to (3.10) of Ref. 0. 
First, when calculating the coefficients of high tempera- 
ture expansion one gets an amplification factor of order 
|r(7 + iu>)\ ~ 5 x 10 10 . Second, while calculating the 
extrapolated slope, one gets another amplification by a 
factor w 3 ~5x 10 3 . Putting these two factor together, 



We are now ready to use the dat a to determine some 
of the unknown quantities in Eq. (5.1). In the follow- 
ing, we will study these corrections for the Ising and the 
LG cases separately. The reason for doing this is that in 
the Ising case, the ratio A\/Aq = —0.56 while in the LG 
case, A1/A0 = 0.14. The relative size of the quadratic 
corrections is presumably of order (Ai/A ) 2 and these 
corrections will be more sizable in the Ising case. We 
start with the assumption that there is one next sub- 
leading correction which dominates when we move from 
bin 13 to smaller values of x. In other words, 



XlO" 



A -^il0 ( 



A1Q- 



(5.2) 



in an intermediate x region. In this equation, the four pa- 
rameters 7, A, Aq and A\ are understood as their best 
estimates near criticality obtained in the previous sec- 
tion. Anticipating the results obtained below, the expo- 
nent 4> is roughly of order one. The corrections in bin 8 
are thus of order 10~ 8 which is precisely of the same order 
as the numerical errors if we use double-precision calcu- 
lations. Consequently we had to use t he Ma thematica- 
based method described in the section III B in order to 



get at least 4 significant digits for the corrections. For 
time considerations, we have limited our calculations to 
10 points per bin. 



If Eq. (5.2) is approximately correct, the logarithm of 
the l.h.s. should be approximately linear in some region 
of x. This quantity is displayed in Fig. [7] for the Ising 
model and in Fig. ^ for the LG model. 
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FIG. 7. \og 10 [\ X W^'' x) - A - AM' Ax) \] 
versus x — —logio((3 c — j3) for the Ising measure. 
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FIG. 8. Iogio[|xl0 ( " 
versus x — —logio(f3 c ■ 



1X) - A Q - AilO c_Aie) |] 
- (3) for the LG measure. 



One sees that in both case the graph is approximately 
linear over a large region of x. Using a linear fit in each of 
the bins of these regions we obtain <f> — 0.82 and A = 0.4 
with a — 9 x 10~ 4 in bin 6 for the Ising model and 
6 = 1.01 and A = 0.4 with a = 6 x 10~ 4 in bin 4 for the 
LG model. Other bins have larger values of a and values 
of which change by a few percent while moving from 
bin to bin. 

In the Ising case, we have w (^) 2 and 4> ~ 2A 
and we interpret this as a second order (or quadratic) 
effect associated with the first irrelevant corrections . In 
other words, A ~ A 2 if we follow the notations of Eq.(5.1) 



. In the LG case (^) — 0.02 is very small and the dom- 
inant effect in the linear region is the analytic correction 
(4> = 1) behavior. In other words, A ~ A a used if we 



follow again the notations of Eq. ( p.l| ) . 

The departure from linearity occurs in its most ex- 
treme way as dips located near x — 2 in the Ising case 
and x = 10 in the LG case. These dips signal the exis- 
tence of effects of opposite signs. A plausible interpre- 
tation of the location of these dips is that they occur 
at values of x where the lO~ 2Ax prevail over the 10~ x 
analytical corrections. A detailed analysis confirms this 
view for the Ising model, which allows us to neglect the 
analytical corrections in bin 13. For the LG model, two 
effects compete in bin 10 which is dangerously close to 
bin 13 where the parameters are fined-tuned (see below) 
and we were unable to get a clear linear behavior after 
one more subtraction. Our most plausible explanation is 
the following for the LG model. Near x — 10, we have 
A a 10- X ~ -A 2 10- 2Ax which implies A 2 ~ 0.016. With 
this rough estimates \^\ « (4 1 ) 2 which is consistent 
with a second order effect. So if this interpretation is 
correct the quadratic effects are about twice the size of 
the analytic corrections in bin 13 for the LG model. 

In summary, we will use the assumption that in bin 
13, the corrections are mostly second order effects and 
we will neglect the analytical corrections. This assump- 
tion is well-obeyed in the Ising case but is just an order 
magnitude estimate in the LG case. 



C. Accuracy of the Previous Estimate 

We are now in position to estimate the effects of the 
next subleading corrections in the calculation of the crit- 
ical exponents reported in the previous section. First of 
all, we notice that by extrapolating the dominant linear 
behavior described in the previous subsection to bin 13, 
we obtain effects smaller than 10 -11 in the Ising case 
and smaller than 10 -13 in the LG case. This justifies 
treating them as small perturbations in bin 13. When 
we fitted the data in bin 13 without taking these small 
effects into account, we made small adjustments in the 
fitted parameters which allowed us to fit the data with a 
precision comparable to the numerical precision. In order 
to get a rough estimate of how much the next subleading 
corrections led us to misestimate the exponents, we can 
linearize the next subleading corrections about x = 13. 
We obtain a change of the "apparent" slope: 



l An 



-10 



-130, 



(5.3) 



The order of magnitude of the corresponding errors on A 
can be estimated by equating the t erm linear in £7 with 
the term linear in 8 A in Eq. ( 4.1C ). This yields: 



\SA\ 



I^10- 13A <57| ■ 



(5.4) 



We insist that this is only an order of magnitude es- 
timate. Plugging numerical values, we obtain \5~f\ = 
3 x 10~ 12 and |<5A| = 2 x 10~ 6 in the Ising case, and, 
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|<?7| = 2 x 10" 13 and |5A| = 6 x 10~ 7 in the LG case. 
In the Ising case, the estimated errors are slightly larger 
than the discrepancies between the values obtained with 
the two measures. In the LG case, the estimated errors 
are slightly smaller. However, larger uncertainties in the 
error estimates appeared in the LG case. If we use the 
largest estimates for the errors, our final result for the 
first method is: 



U{k) = 1. - 0.35871134988fc 2 + 0.0535372882fc 4 



7 = 1.299140730159 ± 3 x 10" 
A = 0.4259485 ±2 x 10~ 6 . 



(5.5) 
(5.6) 



VI. THE EIGENVALUES OF THE LINEARIZED 
RG TRANSFORMATION 

As explained above, the easiest way to calculate the 
critical exponents consists in linearizing the RG trans- 
formation near a fixed point R*(k) specified by the coef- 
ficients a* i. This can be done as follows. First we express 
the coefficients after n iterations in terms of small varia- 
tions about the fixed point: 

a,u = a* + Sa-n.i ■ (6.1) 
At the next iteration, we obtain the linear variations 



Sdn+1,1 = Mi tin Sa n . 



(6.2) 



m— 1 



The l max x l ma x matrix appearing in this equation is 

da n +i,i 



M, 



l.m 



da r , 



(6.3) 



evaluated at the fixed point. 

Approximate fixed points can be found by approaching 
(3 C from below and iterating until the ratio On+i i/a n ,i 
takes a value which is as close as possible to 1. This 
procedure is described in Ref. [Q. The approximated 
fixed points obtained with this procedure depend on (3 C . 
Using their explicit form which we denote R*(k,/3 C ), we 
obtained a universal function U(k) by absorbing (3 into 
k. More explicitly, we found that 



(6.4) 



is in very good approximation independent of the model 
considered. This function is related to a fixed point f(s 2 ) 
constructed in Ref. M by the relation 



(6.6) 



It is not known if there is only one non-trivial fixed 
point for Dyson's model. Using the parametrization of 
Eq. (pTl|), we have considered Q the 12 cases obtained 
by choosing among the following possibilities: m 2 = ±1 
(single or double- well potentials), p — 2, 3 or 4 (coupling- 
constants of positive, zero and negative dimensions when 
the cut-off is restored) and g — 10 or 0.1 (moderately 
large and small couplings). All approximate fixed points 
we ha ve c onstructed give a function U(k) very close to 
Eq. (|6.6|) . The closeness can be characterized by the 
p- norms introduced in [Q. For p = 2 and I < 42 we 
found that the error Sui on the Z-th coefficients of the ap- 
proximate U(k) with respect to the accurate expression 
obtained from Ref. were bounded by \Sui\ < — 
for calculations using double precision. In other words, 
the function U(k) seems to be independent of the general 
shape of the potential, the strength of the interactions 
and whether or not the model is perturbatively renor- 
malizable. 



Using these approximate fixed points, we were able to 
obtain 7 and A with 7 decimal points. In the following, 
we will use directly the more precise function U(k) con- 
structed by Koch and Wittwer [Q. We retained 16 sig- 
nificant digits for the coefficients appearing in Eq. (6.6) 



and used values of l ma x up to 65. We then calculated 



the eigenvalues of the matrix given in Eq. (6.3) with 
two different methods. The first was using "blindly" the 
instruction Eigenvalues in Mathematica. The second 
consisted in using the eigenvalue routine LAP AC K JTJ] 
for which we were able to vary the control parameters 
of the program. The two methods gave identical results 
with 14 decimal points for the first two eigenvalues. The 
first six eigenvalues are given below. 



n 




1 


1.42717247817759 


2 


0.859411649182006 


3 


0.479637305387532 


4 


0.255127961414034 


5 


0.131035246260843 


6 


0.0654884931298533 



f/(fc) (X /((^-)fc 2 ) . 



(6.5) 



The Taylor coefficients of / can be found in the file 
approx . t m 0. Normalizing Eq.([3~J) with {7(0) = 1, 
we obtain 



In order to get an idea regarding the asymptotic behav- 
ior of the eigenvalues, a larger set of values is displayed 
in Fig. |. 



13 



10 



30 



completely calculable. If the procedure can be extended 
to the hierarchical model, then we could almost consider 
the model as solvable. 

Assuming that the non-universal quantities can be cal- 
culated in a reasonably simple way for all the renormal- 
ized quantities, we would be in position to decide if the 
introduction of the bare parameters can be replaced by a 
choice of non-universal quantities appearing in the scal- 
ing laws. If we knew the range of these non- universal 
quantities and their mutual dependence, we could just 
input an "independent set" and obtain directly all the 
scaling laws. In particular, in D — 4, this procedure 
would yield triviality bounds. This alternate way of us- 
ing input parameters in field theory is now being inves- 
tigated. 



FIG. 9. logio[A n ] versus n. 



It is clear from the figure that X n falls faster than ex- 
ponentially with n. This property is important when the 
non-linear effects are calculated. 

Using the first two eigenvalues and the relationship be- 
tween the eigenvalues and the exponents reviewed in the 
previous section, we obtain the values 



7 = 1.2991407301586 ± 10~ 13 
A = 0.4259468589881 ± 10" 13 



(6.7) 
(6.8) 



The (conservative) estimation of the errors is based on er- 
rors of order 10~ 14 on the eigenvalues and the fact that 
the derivatives of the exponents with respect to the eigen- 
values yields factors less than 4. 
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VII. CONCLUSIONS AND PERSPECTIVES 



In conclusion, we have calculated the exponents 7 and 
A with a accuracy significantly better than in Ref. |Q. 
The three independent calculations performed here agree 
on the following value for the leading exponent: 



7 = 1.299140730159± 10" 



(7.1) 



Our results show the excellent agreement between the 
methods developed in Ref. Jn| and an expansion about 
the fixed point of Ref. p|. As far as the calculation of 
the exponents are concerned, the linearization procedure 
is much simpler and more accurate. 

It is important to know if the non-universal quanti- 
ties Aq, Ai, A 2 , A a , . . . could also be calculated by a 
using an expansion about the fixed point which involves 
non-linear terms. We have addressed this question in a 
simplified model, namely the recursion relation for the 
susceptibility 



(7.2) 



A detailed analysis |Tl[] , shows that in this model, the un- 
known quantities appearing in the scaling law for \ are 
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